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Determine f from d intensions 
have i £ hicknesses 
Index of air 

Angle of incidence , degrees 

Incidence angle, radians 

Pel- erw .7 ne the substrate index 

Petermi ne t he FI LM index 
N will always be odd 
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COUFLEDWA VE WL ; TT ; PD 
n Set ORDERS ~ the number 
WA VELENGTH<-WL 
l~< LAYERl ; 2 ] ? OR A T T N O'FFR TOP 
d< LAYERl ; 3 ] 

/)0 4-1 

THETA <-TH 

THET--THETAy-ol ; 18 0 
ns-STTNPEX WA VELENCTH 
n< Op 0 

Fl LHINPEX" LAYERl ; 1 ] 
A'<1 +0RPERG*2 
hH l N ) - 1 
.!>/]-( ( iV-1 ) 5 2 ) 

r\ 

1 <■ ,DP<( N . N )pO 
TTHN.N )p\N 

IT ( 0= , ( 7T -W ) ) / \ N 1 2 I' .1. 
7' C pf )pc(p/)/) ) P 7 
7 7/ <•■ 1 ^ 7 

M.)<o2 i WAVELENGTH 

!<:< i. <■ kO * ( ;j0 v - 1 o 77/ F 7' ) v ?,'/! VE LENGTH f (7i/1 TINGFER1 0D 
klziH ( -( 7T<0 )*2 )+ 1 )x( 7 T<( ( /Y0 + 2 )*( ;;0 * 2 ) ) ■■ ( kxi.*2 ) )* . 5 
!' 7:?:.7>( ( (/<<">* 2 ) ■? ) ) ( /<x i * 2 ) ) '• . 5 n Absorbing 

TM : 

B< ((K+ .*"(EE'®"F<- FERMI TTI VI TY ) )+ . x"A'< M V EN UMBER ) - J 
"n >IF 

EIGENSTUFF E+.y"P n rw eiaenspace calculations 

V<(EE+ .*"W)+ .y"Q n PRODUCT MATRIX FOR TM (EE IS E"E) 

X*Ix" + -k<)y"Qy-d n d SCALAR OR VECTOR WITH LENGTH OF f 

PELT AH ( 2*/V ) , 1 )p( i--0 ) , ( ( ?.oTHET)*Q<n * nQ ) * i~ 0 n 



n 7 .i.s ident ity matrix 



substrate ( Si ) 



F07? 2W 



f ( ( n0*2 )*/<0 ) 

5 ( ( /].S* 2 ) x/K) ) 



Zl+-( l:>7)x(/V,/V)p./<l" 
Z2< ( 1=7 )y(N.N )p!<7z 
W1 <-7 7, , [ 1] ■• 071 *Z1 
FO( 1 ^ 7 ) . I. I ] n 7 1 y.Z2 
FANPG"<b\ p f 

R< N t , ( - DELTA )R)( Ml . - F<v > 
n Pif fraction efficiency for TM 

PERTMH THETAOUT-JH )/( PERTH- 0 )/PERTH<( R*+R )*9o( klzii kQ*n0*2oTHET) 
n PERTH H PERTH* O ) / PERTH '\ R* + R) y- 9 o( .V i n / s W) x ;i0 x 2 oTHET ) 
n 

PF/?7"F' 0 
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A H A'+ • X "A' ) • "'v' 
F/r;FA'.777/7"F 

y- 7," Q 
X ^ 1 y " * - kO v 



7F e iaenspace calcu lat ions 
FROPUCT MATRIX FOR TE 
d SCALAR OR VECTOR WITH LENGTH OF f 

n FOR TE 
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• Q*d r\ 
PELTAH(7yN ) , 1 )p( 0 ) , ( ( ?o77/F7')x07 1.^)0 )x .7>0 
yi*( 1 J )y( N, N )r klr.ii i<n 
Y2< ( 13J )y( N ,N )pk?zi = h 0 
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n Diffraction efficiency for TE 
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E I GEN STU FF EI. 

Z< EI GEN" EI n The funct ion E ' ''-FN is an IBM program product 

W< ( ( of ) p c 1 0)1* Z n cannot l.-e sIikwii here. 

QQ<( ( pf )pc( ( -N ) . 0 ) )4 "S 
<?<-0p0 

<?<-<?x" J 

EIGENVALUE QQ 
Q*-Q.c{N ,N)pQQ* . r > 

FANDG L;XA ;XL;WL: VI 

XL'-L^X 

WL'-L^W 

VL V 

a b <-■ ( i( ( -ft'/., ) . i; 1 1 vr, > . fv; )+.* ( r.7,+ . x,v/, ) . i 

A<(N.N)pAB 

FG < ( . * 77, + X/! ) , i; 1 1 17, + . * TI-XA+-XI* . yA 



VL + . y-XL 



'J 



FILMINPEX FJ LM ; C'1 ; C2 ; €3 ; I 

1 + ( 7. 0 = + / " ( ( r. :fl hM ) ~ CA UCH Y f. ; 1 1 ) ) / i 1 1 p '-A nci-IY 
G\ < CAUCHY[. I ; 7. ] 
C2<-GAUCHY[.I ;3] 
C3<-CAUCHYlli'\ 3 

n^n.CH HC7i( WA VELENGTH y 10)*? ) + C3 s ( .V/1 VELENGTHy-XO )*'* 

E< 'FERN I TT I VITY 
E< OpO 

FFRMFRIME" i o/" 

FERMFRIME M 
FF<-( N , N )p/?+ 1 
7 I < K)FF 

E£ ( ( nUMJ*? ) ■-( nO*? ) )x( 3 o( o 1 x ( J J - FF ) x ZL Af !l ) ) sol y-II-FP 
EEi ( 0 = . ( II -FF) )/ \ A'*? 1 < ( ( n [;/•/;! - ? WTM I ) + ( nO * ? ) x ( 1 - f [ Ml ) 
F> F . r.( o J I )pEE 

K'WAV EN UMBER 
!<<(N ,N)pkxi.-;kO 
!<<■( rK ) y" I 

ns.'-G T. INDEX h'A VE LENGTH : .7 NPEX ; A % ks 
n Determine t h<* complex refract ixe index from ? 1 0 f"o 825 nm. 
7 NPEX<- "! + ?♦( h'A VELENGTHZG II ; 1 J )/i .1 tp,'..7 

ns'Gl[.INPFX\ 1 1 ; ? 1 + ( A ' { h'A V E LEN GTH - G .! I. 1 NPEX\[ .1 ') ; 1. ] ) f^/Gl [ INDEX ; 1 j" ) 
[.INDEX; ?] 
.V s <■• N J 1. 1 N P E X I. t .1 ; 3 J + A y - / G U.I N D E X ; 3 '] 
ns'-ns - 0J \ y-ks 



The function COUFLEPW AVE is called by 
COUFLEPW AVE WL 

where WL is a required argument ; .its value being the wavelength at which 
to evaluate the theoretical profile. COUFLEPW AVE , as configured above, 
is set up to compute TM diffraction. To change to TE , remove the comment 
symbol from line 2 6 & 63 and add a comment to line 62. COUPLEDWAVE_also 
requires several other vari.al.il.es to be defined in tlte workspace: 
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